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Abstract 



We analyze electron transport through a quantum shuttle for the applied volt- 
age below the instability threshold. We obtain current-voltage characteristics 
of this system and show that at low temperature they exhibit pronounced 
steps. The temperature dependence of the current is calculated in the range 
from 2K to 300K and it demonstrates a wide variety of behavior - from 1/T 
decreasing to an exponential growth - depending on how deep the shuttle is 
in quantum regime. The results obtained are compared to experimental data 
on electron transport through long molecules. 
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Recent achievements in nanotechnology facilitate development of a new generation of 
nanodevices incorporating mechanical degrees of freedom, nanoelectromechanical systems 
(NEMS) [1-3]. In this, one of the most promising device structure is a single-electron shuttle 
[4] carrying an electron flow between two contacts. Even classical aspects of this model were 
shown to be important for a description of electron transport through a nanomechanical 
bell [5] and a C eo single-molecule transistor [6]. Theoretical analysis of such a system in 
the classical framework was done in Refs. [7-9]. However, many properties of NEMS can 
only be understood on a fully quantum mechanical basis. Quantum regime of the shuttle 
having mass m and characteristic frequency uo$ occurs when an amplitude of zero-point 
mechanical fluctuations, Ax zp = ^Jh/2muj, is of order of a tunneling length A involved in 
dependence of the tunnel matrix elements on the distance x between the shuttle and the 
leads, Tk a (x) = Tk a exp(x/ \ a ) (a = L,R;Xl = —\,\r = A). Effects of quantization of 
the mechanical motion were examined in Refs. [10-12] and the occurrence of the shuttle 
instability was predicted in Refs. [13,14] for the zero temperature case. 

In the present paper we examine the electron transport through the shuttle system at 
low applied bias voltage (below the instability threshold) with the special stress on the 
temperature effects. To accomplish this, we perform a combined theoretical analysis of 
mechanical motion and electron transport in the model of the quantum shuttle carrying a 
single electron level (strong Coulomb blockade regime). We show that at low temperature 
current- volt age characteristics exhibit pronounced steps. The first step takes place when the 
applied voltage compensates the initial separation between the single electron energy level in 
the shuttle. The second step occurs when the applied voltage facilitates electron tunneling 
from the left lead to the shuttle with absorption of virtual "quantum" of mechanical motion 
(phonon) by electron with consequent emission of this "phonon" during electron tunneling 
from the shuttle to the right lead. Moreover, we demonstrate that there is instability of 
electron transport through the shuttle, when the applied voltage is enough to produce this 
"phonon" absorbed by the shuttle system. 

We show that the character of the current through the system as a function of temper- 
ature is strongly dependent on quantum parameter 

h 

2mdJ A 2 

which defines a relative level of zero-point mechanical fluctuations. In particular, we demon- 
strate that systems with relatively small zero-point mechanical fluctuations (uq < 0.1) ex- 
hibit weak temperature dependence of the current- voltage characteristics I(V, T), whereas 
for the quantum shuttle with vq > 0.4 we obtain the current as the exponential function of 
temperature. It should be emphasized that the model of the quantum shuttle can be of prime 
importance for theoretical explanation of some aspects of electron transport through single 
molecules and self-assembled monolayers [15]. An organic molecule is coupled to the leads 
with elastic links [16] and may oscillate as a whole. These oscillations are reflected in the 
temperature dependence of the current- voltage characteristics of the system [17]. It should 
be noted that experimental current-voltage characterizations of organic monolayers demon- 
strate a wide variety of temperature behaviour - from a weak temperature dependence for 
molecules C8, C12, and C16 sandwiched between Au contacts [18] to an exponential temper- 
ature dependence of a conductance and current that has been observed by D. Stewart et al 
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[19] in experiments with a Langmuir-Blodget monolayer of eicosanoic acid (C20) connected 
to Pt electrodes. The latter results have not yet explained because only no temperature 
dependence (in the case of direct tunneling) or activation dependence (in the case of hop- 
ping transport) were expected. We suggest that the mechanical molecular motion may be 
the reason for such exponential dependence. The difference in measured temperature de- 
pendences probably reflects the difference in connections of the molecules to electrodes and, 
consequently, the different parameters u for these systems. 

To examine the electron transport through the quantum shuttle, we start from the Hamil- 
tonian having the form (h — 1, — 1) 

2 2 2 

H = f- + — ^— + E a + a + 
2m 2 

E E ka 4 a c ka - £(T te c+ a + T* ka a+c ka )e x ' x « . (1) 

ka k 

Here a — L, R, X L — —A, \r = A, and V is a voltage applied to the leads [20]. The equations 
of motion, derived from the Hamiltonian, Eq.(l), are given by 



id = E a - ^ T ka c ka 



x/\ a 
ka 

- fP r J 1 np x/\ a 

ou ka J -' k a'^ k a ± ka u & : 

x + u%x = E -^-(T fcQ c+ a + T* ka a + c ka )e x ' x ". (2) 

ka mA <* 

It should be noted that we take into account nonlinearities associated with a nonlinear 
exponential dependence of the tunneling elements on the distance between the shuttle and 
the leads. 

The electric current through the shuttle is defined as I — I L — —Ir, where 

^ = 4E(4A) = -teJ2TUa + e^c ka ) + h.c. (3) 
dt I I 

With the expression for the electron amplitude in the lead derived from Eq.(2), as 

Cka(t) = c£(f) + iT ka J dhe-^-^ait^^Oit-h), (4) 
and using the formula 

<«W (t)/A -<ffi(0> = 

-iTka\ dti/«(^*«)e- iBh ^*-* 1 ^(t-ti)([a + (0e*W/ A «,o(t 1 ) e *(* 1 )/ A -] + >, (5) 

resulted from the Langreth theorem, we obtain the electric current, as given by 

I a = -eY,\T k a\ 2 I dhUE^e-^-^eit - h){[a + {t)e*W^ , a(t 1 )e^ (il)/A "] + ) + 
k J 

El T ^| 2 / ^ie-^ (i -* 1 ^(t-t 1 )(a+(t)a(t 1 ))(e a; W/ A ",e a;( * l)/A ") + h.c. (6) 



Electrons in the a— lead are assumed to be described by the Fermi distribution f a (E) = 
f(E—/j, a ) = exp [ E " r ^ a ^j + 1 . For averaged values of electron correlators we have approx- 
imate expressions (a+(t)a(ti)) = Ne iEo(t - h) , (a(ti)a + (t)) = (1 -iV)e i - Bo( *~* l) with N = (a+a) 
being a steady-state electron population of the shuttle that can be found from the condition 
Ir + Il = 0. Using previously derived formulas (see Appendix in Ref. [17]) in the wide-band 
limit we obtain the expression for the current, as 



+oo +oo 



I a = e Uc eT a J2 X] J m (vo)Ii(v c ){N[l - f a (E - rruv - Iuj )\ - 

m=—oo l=— oo 

(l-N)f a (E + muj + luj )}. (7) 

Here, J m {va) and Ii(v c ) are usual and modified Bessel functions of order m and /, respec- 
tively, T a = 2nJ2k \Tk a \ 2 S(uj — E ka ), and u c = (x 2 )/\ 2 is a relative dispersion of mechanical 
oscillations for the nonequilibrium case. In the following, we consider symmetric coupling 
between the shuttle and the leads, Tl = Tr = T. A steady-state value of the electron 
population in the shuttle is determined by the relation 



with 



N = C/D, (8) 



C = Jm(vo)Ii(v c )[fL(E + muj + + f R (E + muj + lu )] 

ml 

and 

D = ^2 ^m( z/ o)^( z/ c)[2 - fh{E G - muo - lu ) + f L (E + mw + lu ) - 



ml 



f R (E - mujQ - Iujq) + f R (E + muj + lu Q )]. (9) 



In the case of relatively small zero-point mechanical fluctuations, u <^ 1, J2i Ii(v c ) = e^ c , we 
have the simple formulas for the shuttle population iV and for the current /, as 

N = e c }^li{isc) , 

i z 

I = h = ere- E 7|M /^. + M-A(g, + M . (10) 

I Z 

It should be noted that effective temperature of the biased nanomechanical system, which 
is determined by the dispersion of mechanical fluctuations, (x 2 ), can be different from the 
equilibrium temperature T [17], so that we cannot neglect nonlinearities caused by v c even 
when !/ <l. 

To calculate a relative dispersion of mechanical fluctuations, v c = {x 2 )/X 2 , the equation 
(2) for the shuttle position can be rewritten in the form of quantum Langevin equation 

x + 7± + LU 2 x = £, (11) 

with a damping rate 



A 



7 = u Te^B(cu ), 



(12) 



and a fluctuation source, 

e = E(^/^){« + ^ /Aq ,cS} + &.C., (is) 

which is characterized by the spectrum K{uj) = u e Uc T(hujQ/m)A(uo). Here, we introduce 
the following dimensionless combinations of the Fermi distributions, A{uj) and B(u), as 

a ml 

{N[2 - f a (E - mu - lu -uj)- f a (E - mu - lu + uj)] + 
(1 - N)[f a (E + muo + luj + uj) + f a (E + muj + - ^)]} 

and 

B ( u ) = EE J '»W / 'K) x 

a ml 

{N[f a (E - mujQ - lu) - ui) - f a (E - muiQ - Iuj + to)} + 
(1 - N)[f a (E + muo + luo Q - uo) - f a (E + mLU + lu + uj)]}. (14) 

The relative level of mechanical fluctuations, is c , can be found from the self-consistent equa- 
tion: 

Vc = Uo bWy (15) 

supplemented by the equation (8) for the steady-state electron population N of the shuttle. 

We have solved self-consistent set of equations, Eqs. (8) and (15), and substituted the 
obtained values of N and v c into Eq. (7) for electric current. Corresponding current-voltage 
characteristics are shown in Figure 1 for various temperatures, energetic separation between 
the electron level in the shuttle, Eq and equilibrium chemical potential of the leads, /x, being 
0.2uj , u)q/2ii = ITHz, and for quantum parameters u = 0.7 (Fig.l(a)) and v> Q = 0.07 
(Fig. 1(b)). It is evident from these Figures that there are pronounced steps in the low- 
temperature current- voltage characteristics at the bias voltages (i) eV/2 = E — /i and (ii) 
eV/2 = uj — Eq + fx. The first one occurs when the chemical potential of the left lead 
passes through the energetic level of the shuttle, while the second one corresponds to the 
passing of the chemical potential of the right lead through the virtual level with energy 
E — uq. In the latter case, tunneling of an electron of the left lead having energy E — uj q 
to the shuttle level with energy Eq is accompanied by absorption of virtual "quantum" of 
mechanical motion (phonon) with further tunneling of this electron to the state of the right 
lead having energy Eq — uoq accompanied by emission of this "phonon". This second step is 
much more pronounced for the case of larger v . When the chemical potential of the left lead 
passes through the virtual level with energy E + u Q (at the voltage eV/2 = u + E — fx) , 
electron tunneling from the left lead to the shuttle can be accompanied by absorption of 
such a phonon by the shuttle, and the system becomes unstable. To illustrate this, we 
present in the Inset of Fig. 1(a) the current- voltage characteristics (in the logarithmic scale) 
with various values of E — fx. It should be noted that the condition for this instability is 
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identical with that found in Ref. [14]. Finally, one can see in Fig. 1(a) that current- voltage 
characteristics become more smoother and steps disappear with temperature increasing. 

The temperature dependence of the current through the quantum shuttle system is of 
special interest. It is shown in Figure 2 in the logarithmic scale for Eq — /j, — 0.2uio, eV = 
0Au> , ujq/271 = ITHz, and various values of quantum parameter v . It is evident from this 
Figure that initial decrease of the current with temperature increasing is followed by the 
strong exponential growth at larger values of u . The exponential temperature dependence of 
the current- voltage characteristics has been also predicted for the model of a tunnel junction 
having its matrix elements modulated by the vibrational motion [17,21]. In this model [22], 
the position of the oscillator affects directly a probability of electron tunneling between leads 
and the temperature dependence involved in distribution functions of the leads is integrated 
out. Correspondingly, the exponential temperature dependence, caused by fluctuations of 
an oscillator displacement, occurs for any levels of zero-point mechanical fluctuations, and, 
in particular, for z/ 1. The present study reveals that an inclusion of a single-electron 
state into the vibrating molecule leads to an additional factor 1/T in the formulas for current 
and conductance (caused by the temperature dependencies of the distribution function of 
electrons in the leads) that smooth out the exponential T-dependence in the case of small 
uq. Indeed, for high temperatures, (T ^> \E — fi a \, one can see that v c = i>o(2T/huJo), so that 
the current, Eq.(7) is described by the formula 

eV 

I = eTexp(4Tv /tw ) — , (16) 

where a temperature dependence reveals itself not only in the exponent e /kTu °^ hul0 as for the 
tunnel junction coupled to the mechanical oscillator [17,21], but also in the factor 1/T. An 
outcome of a competition between these two factors depends crucially on the value of the 
relative level of zero-point mechanical fluctuations u Q . At small u the factor 1/T dominates, 
whereas at larger z/ the current-temperature curve demonstrates a very pronounced expo- 
nential temperature dependence. However, u depends on the inverse tunneling length of 
system squared and even small variations of the tunneling length can bring about completely 
different behavior. Accordingly, the experimental data of Refs. [18,19] might correspond to 
slightly different conditions of connections between molecules and leads and still exhibit very 
different temperature dependence of the electron current. 

In conclusion, we have analyzed electron transport through a quantum shuttle having 
single electron energetic level. We have shown that at low temperature there are pronounced 
steps in current- voltage characteristics corresponding to direct tunneling through the shuttle 
and to tunneling accompanied by consequent absorption and emission of virtual " quantum" 
of mechanical motion (phonon) by electron. We have found that when the applied voltage 
facilitates the emission of "real phonon" during electron tunneling from the left lead to the 
shuttle, the instability of electron transport occurs. The temperature dependence of electric 
current has also been determined in the range from 2K to 300K. Mechanical motion of the 
shuttle leads to the exponential temperature dependence while the temperature dependence 
of the electron distribution function in the leads brings about the factor 1/T at high temper- 
ature. The competition of these two factors gives rise to wide variety of current-temperature 
curves depending on quantum parameter z/ , such as decreasing of current with temperature 
increasing (at small z/ ), a weak temperature dependence (at the intermediate values), and 



a strong exponential growth (at large u ). This variety was previously demonstrated in the 
experiments on electron transport through long molecules [18,19]. 
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FIGURES 



FIG. 1. Current-voltage characteristics of quantum shuttle system at Eq — fi = 0.2ujo and at 
various temperatures; (a) for quantum parameter vq = 0.7, (b) for quantum parameter vq = 0.07. 
Inset: Instability of electron transport at higher voltages for various separations of the electron 
energetic level in the shuttle and equilibrium chemical potential of the leads. 

FIG. 2. Temperature dependence of the electric current through the shuttle for various values 
of quantum parameter vq. 
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